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Weak Magnetic Order in High-Tc Superconductors Produced by 
Spontaneous Josephson Currents 



. -E 

rn ■ 
(N 

<: 

oo 



o 
o 

u 

§-: 



o 



V. L DE Mello^'^ and David Mockli^ 

Department of Physics, Simon Eraser University - Burnaby, British Columbia, Canada V5A 1S6 
Institute de Fisica, Universidade Federal Fluminense - Niteroi, RJ 24210-340, Brazil 



PACS 74 . 20 . Mn - Nonconventional mechanisms 
PACS 74. 81. Fa - Josephson junction arrays and wire networks 
PACS 74 . 62 . En - Efltects of disorder 

PACS 74.81.-g - Inhomogeneous superconductors and superconducting systems, including elec- 
tronic inhomogeneities 

Abstract - We develop a model for high-Tc superconductors based on an electronic phase sepa- 
ration where low and high density domains are formed. At low temperatures this system may act 
as a granular superconductor forming an array of Josephson junctions. Cuprates are also known 
to have low superfluid densities and strong correlation effects. Both characteristics activate a neg- 
ative Josephson coupling due to frustration that leads to spontaneous currents responsible for the 
weak ferromagnetic order. This original approach reproduces the observed onset of spontaneous 
magnetic signal and its dependence on the doping level. 
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Introduction. — Understanding of the normal state 
'pseudogap" phase of higli-tcmperature superconductors 
Eemains one of the major puzzles of condensed matter 
jihysics. The reason abides in unusual properties and par- 
tial gapping of the low energy density of states below a 
dertain temperature T* [iHS] , taken as the boundary of the 
pseudogap phase. The lack of an accepted theory hints the 
ilanoscale complexity and the intrinsically inhomogeneous 
electronic structures [HIS], which might vary according to 
the specific family of cuprates. 

• Among the many anomalies of the pseudogap phase 
[T}{3], the weak ferromagnetic signal on YBa2Cu306+a; by 
zero-field muon spin relaxation [5] and by polar Kerr ef- 
fect [7] (PKE) requires an explanation. Both measure a 
hole concentration (p) dependent signal that starts jointly 
with T* (p) in agreement with many different experiments 
[IHa]. The values of T*{p) are well above Tc{p) in the 
underdoped region and drops rapidly with increasing p, 
becoming comparable with Tc near the optimally doped 
concentration p — 0.16. 

In this letter we address this long standing problem by 
an entirely original approach. We assume an electronic 
phase separation with the charges segregated in domains 
(grains) of low and high densities separated by a potential 
barrier [5HTU] . This charge segregation reduces the kinetic 
energies enhancing the possibility of superconducting am- 



plitude formation in the isolated grains as a granular su- 
perconductor. This local superconductivity in domains is 
possibly similar to the case in which Bi-clusters of 2.5 to 
40 nm becoming superconducting, in contrast with the 
non superconducting bulk Bi |11| . The free energy barri- 
ers between the domains produces charge confinement, as 
in the Bi-cluster case, and acts as weak links. We model 
this system as an array of Josephson junctions. Low local 
charge densities yield strong phase fluctuations [T^ that 
combined with strong correlation effects gives rise to frus- 
tration [13| . Our breaking new idea proposes that these 
frustration effects are the origin of the weak ferromagnetic 
signal in YBCO [HIIT], which is a generalization to single 
crystals of the method used to explain the paramagnetic 
response of granular cuprates [T^ . 

Many reviews j4j[5] have discussed the existence of elec- 
tronic inhomogeneities by different experiments on several 
samples. The disorder of the charge in cuprates can be 
treated by the Cahn-Hilliard differential equation that de- 
scribes the formation of patches or grains. The charges 
trapped in these grains lose kinetic energy, form bound 
states and Cooper pairs, give origin to the pseudogap and 
superconducting phases [Hllin]- With this, we illustrate a 
typical Ginzburg-Landau free energy simulation [5] on a 
105x105 square lattice in Fig.([T]). 

Earlier experiments on granular cuprate superconduc- 
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Fig. 1: (color online) Typical simulation of the Ginzburg- 
Landau free energy potential VoLiu, T) on a 105 x 105 square 
lattice. The grain's diameter spans a few nanometers. At low 
temperatures the charges are attracted to the grains, similarly 
to a granular superconductor. We show the d-wave amplitude 
on the top left. 
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Fig. 2: A schematic example of a sequence of intermediate 
states for a Cooper pair tunneling across a wall which contains 
a spin up and leads to a negative coupling along the junction. 
Between steps b and c it is necessary to permute the two elec- 
trons to be in the canonical order of a pair, and this exchange 
is responsible for the negative sign |13) . 



tors showed anomalous magnetic properties attributed to 
frustration due to properties of the Josephson junctions 
[H] . In this case, the frustration effects arise due to the 
crystal's axis mis-orientation of the grains with supercon- 
ducting d-wave order parameter; this gives a negative con- 
tribution to the Josephson-j unction energy. The current 
is negative with a phase shift tt -a tt junction. Sigrist and 
Rice [13] demonstrated that a loop of current with an odd 
number of tt junctions is frustrated, i.e., there is no way 
to minimize the free energy of all junctions in the loop 
and a spontaneous current may arise if the coupling is 
sufBciently strong. 

However, we cannot talk of mis-alignment of the crys- 
tal axis between different grains on a single crystal. In 
such system, negative Josephson coupling is possible due 
to the large fluctuation in the local density and correlation 
effects. This was demonstrated in detail by Spivak and 
Kivelson [T3] using perturbation theory up to 4*''-order 
in a model where correlations effects produce a negative 
Josephson coupling across a junction. For completeness 
we illustrate one of the situations described in their pa- 
per where a Cooper pair tunnels through a barrier which 
contains a spin up as shown in the Fig.©. 

It is important to remark that the physics of d-wave 
Josephson junctions presents many different and interest- 
ing aspects, depending on the type of superconductors, 
order parameter symmetry and geometry of the interface 
[23) . The d-wave pairing order parameter is very sensitive 
to inhomogeneities and interfaces. Quasiparticle scatter- 
ing at interface distorts the order parameter and influences 
the Josephson effect as well as the quasiparticle tunnel- 
ing current |24j . This phenomenon can generate midgap 
states, i.e., surface or interface states with zero energy rel- 
ative to the Fermi energy j25j or Andreev zero-energy. The 
midgap states enhance the Josephson current near T = 



and are observed through the zero bias conductance peak 
P5H27] . These states are robust phenomena and repre- 
sent a crossover from a traditional Josephson junction (0- 
junction) to a 7r-junction, or a phase with broken time 
reversal symmetry. 

On the other hand, the very high temperatures up to 
200K of the signal measured by Xia et. al. [7| shows that 
its origin is not likely to be due to these interface states. 
Furthermore, the intrinsic inhomogeneous nanometer do- 
mains have small and randomly oriented interface, in 
opposition to those specially prepared junctions p5H27j . 
and consequently it is difficult to generate midgap states. 
Thus, here we explore the idea that the origin of the tt- 
junctions and the spontaneous magnetization in YBCO 
is due the very large phase oscillations due to very small 
local density of Cooper pairs [T^ . 

Calculations. — The Cooper pairs tunneling through 
a localized state in the junction between two grains ex- 
change order; this generates the negative sign -a charac- 
teristic of a TT junction |13j . Since the electronic phase 
separation with the granular structure provokes strong 
charge fluctuations, -this- and not the mis-orientation 
of the grains, is the mechanism originating tt junctions 
throughout the system. 

Following the procedure of Sigrist and Rice [H], the 
properties of the whole loop of an odd number of tt junc- 
tions is given at one weak link. Assuming that the current 
/, which flows in the loop, is small compared to the critical 
current of the grains, wc can write |14) 



F(/,A0) 



LP 
2c2 



27rc 



cos (A0 -I- a) . 



(1) 



where A^ is the phase shift across the junction, L de- 
notes the self-inductance of the loop, Ic the critical current 
through the junction, and $0 is the flux quantum hc/2e. 
The additional phase shift a is tt or whether the loop 
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is frustrated or not. Assuming also that the current flows 
through one point, this leads to 



A(/) = 27rn - 27r 
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where n is an integer, $ is the total flux threading the loop 
with three types of contributions; the external field, the 
flux from other loops called ^ex and the current /. With 
this A(j) EqlH becomes [13] 



F{I, $ex) 



LP 
2c2 



27rc 



2^ f^ LI\ 



(3) 



By minimizing F with respect to / for a given ^ex, we 
relate / and ^ex- At zero external field we are interested 
in the case where ^ex ~ because the contributions from 
the other loops should be small. Then, for a = tt, a spon- 
taneous flows around the loop if 7 = 2ttLIc/^oc is larger 
than 1 (Fig. ([3]). This spontaneous current -due to the 
loop's frastration - generates the observed weak magnetic 
signal. 
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Fig. 3: The induced current in a frustrated loop for three 
values of the parameter 7 as function of the external field. 
For the zero field case ($ = 0), it shows that for 7 > 1 a 
spontaneous current may arise while there is not any for 7 < 1. 



that a spontaneous current arises. On the other hand, 7 
is proportional to the critical current Ic, and this depen- 
dence was studied in detail for a d-wave superconductor 
by many authors [iGlfTTj . 

In particular Bruder et al |17j calculated the supcrcur- 
rent tunnel matrix elements in second-order perturbation 
theory for two d-wave superconductors (I and 2) with su- 
perconducting amplitude 

^d.i/2{h T, </)) = A(z, T) cos [2(0 - 01/2)] , (4) 

where is the azimuthal and 4>i/2 is the mismatch an- 
gle which is here zero because the electronic domains are 
always aligned with the crystal axis. The dominant con- 
tribution to the critical current is from "node to node" 
tunneling and the overall behavior is like an s-wave su- 
perconductor [17] . This result is in agreement with the 
calculations of a Josephson junction with two different s- 
wave superconductors |18] . which yields a supercurrent 
proportional to AiA2/(Ai -I- A2). 
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Fig. 4: The temperature evolution of Ad(i) at five differ- 
ent locations of an average p — 0.15 compound. The mean- 
field BdG calculations all vanish at the same temperature 
T*(0.15) ^ 165K (above rc(0.15) = 92K). The low temper- 
ature values between 30-66meV are in the energy range of the 
Bi2212 LDOS gaps measured by STM. 



Thus, as the parameter 7 surpasses unity, a spontaneous 
current flows around frustrated loops and local magnetic 
signals arise. The first spontaneous currents generate a 
magnetic field that infiuences the direction of other cur- 
rents and small local inhomogeneous ferromagnetic signals 
arise throughout the sample, similar to the zero field sus- 
ceptibility of granular superconductors near Tc [15] (T* 
in our case). Muon spin relaxation (/x-SR) [^ and polar 
Kcrr-effect [7] measured such magnetic signal on YBCO 
compounds. 

To compare our proposal with these experimental re- 
sults we need to study the temperature and doping de- 
pendence of the parameter 7, since it controls the sponta- 
neous currents. 7 needs to be larger than unity in order 



Bogoliubov-deGennes calculations demonstrate that the 
values of the local superconducting amplitude also vary lo- 
cally inside a system with local electronic disorder [1"9H2T| . 
In Fig.dl]) we show a typical result of these calculations and 
plot the temperature evolution of the superconducting am- 
plitude Ad{i) at five randomly chosen different locations 
or domains i. 

The disorder of the local gaps implies different values 
of the tunneling matrix between different grains. Conse- 
quently the values of 7 also change at different junctions 
(7(1)). The values of 7(1) increases jointly with A(i, T) as 
the temperature decreases. Therefore the onset of spon- 
taneous current occurs when the largest A{i,T) reaches a 
critical value causing some junctions of a certain loop(s) to 
have ^max > 1 at zero external field, as shown in Fig.Q. 
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Fig. 5: The temperature evolution of the maximum amplitude 
of the d-wave gap A{p,T) of six different values of average 
doping p. The line at ~ 40nieV shows the experimental onset 
of spontaneous magnetization and the intersections yield the 
temperature onset for different compounds. 




Fig. 6: The calculated onset of spontaneous magnetization as 
derived from the intersections in Fig. ([5]) of the A,nax{p, T) with 
the 40meV horizontal line together with the the experimental 
results of /^-SR [6] (open circles) and polar Kerr-effect [7] (black 
squares). The line connecting the theoretical points are just a 
guide to the eyes. 



Results. — It is very difficult to calculate theoreti- 
cally the actual values of 7(1) but we can obtain some 
estimation from the experimental data. For instance, it is 
reasonable that A(i,T) is the major parameter that con- 
trols 7(1) as the hole doping p is varying. The potential 
barrier between the grains is important but it is also con- 
nected with the values of the superconducting amplitude. 
Underdoped systems arc in general more disordered than 
the overdoped ones, then it is possible that the number of 
loops is larger in this limit increasing the strength of the 
ferromagnetic signal. Indeed the PKE signal observed in 
underdoped samples are much stronger than those mea- 
sured near optimally doped [22]. Furthermore, below Tc 
the phase coherence diminish the phase oscillation causing 
the signals measured below Tc to be so small that can be 
almost attributed to oscillation of the data [HIIZ] . 

Since the magnetic signal vanishes near p ~ 0.18, we 
take the onset of A(i, T) that triggers the spontaneous cur- 
rent as the value shown in Fig. ([5]) of Amax{p = 0.18, T = 
0) « 40)meV. In other words, a spontaneous current can 
exist only above A « 40meV. 

With this maximum value taken as the onset to have 7 
larger than 1, we can follow the temperature evolution of 
the maximum gap of each compound with the line drawn 
at Amaxip = 0.18, T = 0) (Fig.®). The results for 5 
doping values are shown in Fig. ([5]) together with the the 
experimental results of /i-SR [Q (open circles) and of polar 
Kerr-effect |3 (black squares). 

Conclusions. — In summary, we used the fact that 
cuprate superconductors have an intrinsic electronic dis- 
ordered state where the charges are segregated in a few 
nanometers high and low density grains separated by thin 
potential walls. As in Bi-clusters [TT] this charge confine- 
ment may be the origin of the local superconducting in- 



teraction. We then calculate the superconducting proper- 
ties by a BdG method using a phcnomenological two-body 
potential proportional to the energy barriers or walls be- 
tween the grains [20] . The various distinct regions are 
coupled forming an array of Josephson junctions that pro- 
mote the resistivity transition at low temperatures. The 
phase-number quantum fluctuations are large enough and 
strong correlation effects promote frustration or negative 
Josephson coupling [13J that leads to spontaneous currents 
and an overall magnetization. 

This approach connecting an electronic disordered state 
forming an array of Josephson junctions at low tempera- 
tures with frustration effects provides an interpretation to 
the intriguing spontaneous ferromagnetic phenomenon on 
cuprates, a long standing problem. This novel calculation 
shows the importance of the charge inhomogencities and 
leads to a pseudogap phase with isolated regions of non- 
vanishing superconducting amplitudes without long range 
order. 
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